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1. Introduction and statement of results 

The solution of a problem arising in integrable systems JUj requires the asymp- 
totics as a —>■ oo, with very small error, of the inverses and determinants of truncated 
Wiener-Hopf operators W a (a) acting on L 2 (0, a), both in the regular case (where 
the Wiener-Hopf operator W(a) on L 2 (R + ) is invertible) and in singular cases. This 
paper treats two cases where a has simple Fisher-Hartwig singularities, one double 
zero or two simple zeros. 

We first state a result for the regular case. Assume a belongs to the Wiener 
algebra of functions of the form a + k with k G L 1 (R), 1 and that 



oo 

<r£)^0, Aarga(0 =0. (1) 

OO 



Then 



W a (a)- L = P a W(a)- L P a - Q a H(1/Z1) H(l/a+) Q a + o(l). (2) 



The notation here is the following. The operators P a is multiplication by X(o, a ) 
or extension by zero from (0, a) to (0, oo), depending on the context. For functions 
defined on (0, a) we define (Q a f)(x) = f(a—x), we use H to denote Hankel operators 
as usual, and for a function v we define v(l;) = v(—£). The term o(l) denotes a family 
of operators whose operator norms tend to zero as a — > oo. The functions a± are 
the Wiener-Hopf factors of a: their product equals a, and a^ 1 resp. crj: 1 extend to 
bounded analytic functions in the lower resp. upper half plane. 

If in addition logo - G L 1 and Jf^, \x\ \k(x)\ 2 dx < oo then the Kac-Achieser formula 
holds: 

det W a {a) ~ G{a) a E{a), (3) 
The Fourier transform we use is 

HO = f-oo &iXi k( - x ) dx. For the Toeplitz case with a weaker 
assumption see Th. 2.14]. For @ itself see footnote [SJ 
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where 

G(a) = exp {i- jT log a(0 df} , £(<r) = exp (log^))' log MO df} ■ 

In the present paper we assume that 

a(0 = (£ 2 "P 2 )r(0, 

where p is real, and r is nonzero, has index zero, and its Wiener-Hopf factors r± are 
of the order as £ — > oo. 

In order to minimize technical details we make a very strong assumption, namely 
that e 5 ^k(x) G L 1 (R) for some 5 > 0. 2 This will result in an exponentially small 
error in the approximation. 

We make a choice of which Wiener-Hopf factor of a incorporates which zero, and 
take the factors to be 

mo = <(e-p)r_(o, mo = -<(e+p)7+(0- 

As functions, we have 

+«-«'• ^-^^ +«+<?), 

<r_(0 i(f-p) MO -hC+P) 

where u± are bounded smooth functions. In fact we think of cr±(£) -1 as distributions 
defined by 

Mp) _ 1 u ^ Mp) _ 1 u ^ 

Observe that if we set e p (x) = e ipx then (0 — i(£ + p))" 1 is the Fourier transform of 
X + e p and (0 + i(£ — is the Fourier transform of X _ e p . (We write X 1 * 1 for Xr±.) 

To state our first result for p ^ we define 

( = H(l/a + )e p , V = H(l/^)e p , 

(here e p is thought of as restricted to R + ) and 

lip 



T-(p) r + (—p) e tap — r + (p) r_(— p) e tap 



2 We shall consistently use S to denote some positive quantity, different for each occurrence. 
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Theorem 1 (p / 0). There is a 5 > such that if A = 0(e 5a ) then W a (a) is 
invertible for large a and 

WaOr)" 1 = P a W(l/a+) W{1/ a.) P a - Q a H(l/al) H(l/a + ) Q a 

+A [e tap r+(p) r_(-p) P a C ® P Q 77 - r+(p) r + (-p) P Q C ® Q«C 
-r_(p) r_(-p) Q a r/ <g> P a r/ + e tap r+(p) r_(-p) Q Q r/ <g> <5aC]- 

This requires lots of explanation. First, the sign = between two operators indicates 
that the difference is an integral operator whose kernel is uniformly 0(e~ 5a ) for some 
5 > 0. 

There are tensor products in the statement. For functions v, w we define v <g> w 
to be the operator taking a function / to v (w, f). (The inner product is the integral 
of the product — no complex conjugate.) In general we have V\ <g> W\ ■ v 2 <8> W2 = 
(wi, v 2 ) Vi <S> W2 and (v <S> w) ■ T — v <S> T'w, where T' is the transpose of T. 

Next, we define W((0 — +p))~ 1 ) to be convolution by X + e p on (0, 00). This is 
defined on any locally integrable function since 

((X + e p )*f)(x) = e^ f X e-^f(y)dy. 

J 

Our assumption on o implies that u + is the Fourier transform of an exponentially 
decaying function, 3 so W{u + ) is defined on any function of at most polynomial growth. 
This extends the domain of W(l/a+) to any function of at most polynomial growth. 
Similarly we define W((0 + i(£ — p))^ 1 ) to be convolution by X~e p . This is defined 
only for functions in L 1 since 

roc 

((X-e p )*f)(x)=e*° / e-*»f(y)dy. 

Thus W(l/a-) acts on L 1 (R+). 

Finally, there appear Hankel operators whose symbols are distributions. The 
symbol of if (1/(7+) is the Fourier transform of a constant times X + e p plus a rapidly 
decreasing function. Thus the operator has a singular part which can be written as 
a constant times e p <g> e p . Similarly for H{l/al). In the statement of the theorem we 

3 Our assumtion implies that <J+{£) extends analytically to a strip around the real line and that 
er+(£ — iS)^ 1 is in the Wiener algebra for sufficicntlly small 6 > 0. Therefore, since (0 — i(£,— iS))^ 1 = 
— (S + it;)^ 1 is in the Wiener algebra, so is — i5). This implies that u+ is the Fourier transform 
of a function with bound 0(e~ Sx ). 
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see the product of two of these, or one of them acting on e p . We interpret these by 
going to the Fourier transforms. For example, from the general identity 



oc 



v , w ) = tt / 



DC 



2tt 

we obtain 

(X + e p , X + e p ) = — j Q + i ^_ p) 0-i(£ + p) ^ = ~ 2ip~> 

since by going into the upper half-plane we pass the pole of the first factor at p + Oi. 
Similarly we find that if((0— z(£ +£>) _1 ) e p = —e p /2ip. (Everything can be comuputed 
by first giving p a small positive imaginary part and then passing to the limit.) 

Although we shall be using distributions throughout they will all be of a simple 
form: Fourier transforms of linear combinations of e± p (combinations of 1 and x when 
p — 0) and exponentially decaying functions, and their restrictions to R . Thus the 
meaning of the operators v — *■ v± in Fourier transform space is clear: take the inverse 
Fourier transform, multiply by X ± , and then take the Fourier transform. 

The statement for p = is simpler since Hankel operators with distribution symbol 
do not appear. We define 




B = a + i 

\ a - ) 

where the factors are chosen so that r_(0) = r + (0). 4 

Theorem 2 (p = 0). The operator W a (a) is invertible for sufficiently large a and 
W a {aY l = P a W(1/<t+) W{l/aJ) P a - Q a H{u^) H(u+) Q a 
- B- 1 [Q a H(ul)l + P a W(l/<r + )l] <8> [Q a H(u+)l+P a W(l/al)l]. (4) 



There are several results in the literature giving first-order asymptotic results for 
the inverse when the symbol has a single singularity or zero. The method of [S] gives 
a first-order result in terms of convergence in operator norm. Rambour and Seghier 
jH] considered the Toeplitz case for a symbol with a zero of arbitrary even order 
and obtained a stronger result in these cases. Ehrhardt |B] obtained results in the 
Toeplitz case for a general class of symbols with several Fisher-Hartwig singularities, 
with convergence in the operator norms of weighted £ 2 spaces, although they do not 
apply to symbols with two real zeros. In this work, too, distributions played an 
important role. None of these serves our present purpose, which is to find formulas 
that hold uniformly throughout (0, a) with very small error. 

4 Without this normalization the second term is to be multiplied by (t_/t+)(0) = — (a-/a+)(0). 
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An earlier and more involved derivation of (j3J than the one we give here is in 
The present proof will use a uniformity in the lemma to Theorem 1 and then 
a passage to the limit. The computation would be tedious if done by hand and we 
resorted to computer computation. For the skeptical reader we give a direct proof in 
§6, along the lines of the proof of Theorem 1. 

For the determinants we assume that in addition logo" 6 L 1 . We use = here to 
indicate that the difference of the two sides is 0(e~ Sa ) for some 5 > 0. Now we define 



G(a) = expj— £^(log<7 + + logcr_)^| 
f 1 f°° 

E(a) = exp j— j^Jloga+Y loga_ d£ 



In the integrals logcx + is interpreted as log (£ + p + Oi) plus a regular function while 
log<r_ is log (£ — p — Oi) plus a regular function. These depend on our choice of factors 
a±. 

For p the result is 



Theorem 3 (p / 0). We have as a — > oo 

det W a (a) = r_(p) r+(-p) - r+(p) r_(-p) e 2iap 
G(a) a = r+(-p)r_(p) 

For p = the result is 
Theorem 4 (p = 0). We have as a — > oo 

det W a (cr 



(5) 



= BE(a). (6) 



This may be rewritten 



det W a (a) =aE{a)+ ./a ± \ + 



G(a) 



<7_ 



Thus it gives the second-order asymptotics with very small error. The first-order 
asymptotics were obtained by Mikaelyan [7j. In the special case er (£) = £ 2 /(l + 
one has 

W a (a )=e- a (l + a/2) 

exactly, as was also shown in Exact formulas for arbitrary rational symbols were 
obtained by Bottcher [2] and used by him to obtain further asymptotic results in 
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these cases. For the Toeplitz analogue the asymptotics are known for symbols with 
any number of zeros [3]. (Or see jH §10.47].) For Wiener-Hopf operators the case of 
two simple zeros was considered in pQ, but the result was obtained only for r real and 
positive. In §5 we show that our formula agrees with the one obtained there. 

Here is how the approximate inverse in Theorem 1 is obtained. Since (J2J) holds 
for regular symbols we consider 

P a Wia- 1 ) WiaZ 1 ) P a - Q a H(l/*1) H(l/a + ) Q a (7) 

a possible first approximation for singular symbols. So we see what happens when we 
multiply it by W a (a) on the left. We find that, with small error, the result is I plus 
an operator of rank two. The form of this operator suggests the ingredients of a rank 
two operator that should have been added to (J7J) to obtain / in the end. A simple 
computation tells us exactly which operator to take, and the rest is straightforward. 
For the determinants we introduce a parameter, use the formula for the logarithmic 
derivative in terms of the inverse, and apply the earlier results. 

2. Proof of Theorem 1 

Lemma. Define 

T a = P a W(l/a+) Wil/aJ) P a - Q a H(l/al) H(l/a + ) Q a 

+A [e iap T+ (p) r_(-p) P a (®P*V- r+(p) T+(-p) P Q (® Q a C 

-T-(P) r-(-p) QaV ®P a V + e iap r + (p) r_(-p) Q a ( ® Q a rj] . (8) 
Then if 6 is small enough we have, uniformly in p, 

W a {a)T a = I + 0(Ae~ Sa ), T a W a (a) = I + 0(Ae~ &a ). 

Proof. We shall use, innumerable times and without comment, the identities 
Win) W a (a 2 ) = a 2 ) - P a H{a 1 ) H(a 2 ) P a - Q a H(ai) H(a 2 ) Q a , 

and 

H{a x a 2 ) = W{a{) H{a 2 ) + H{a{) W(a 2 ), 

and their familiar consequences. (For the Toeplitz analogues of the stated identities 
see H §2.17 and §7.8].) 

Let us multiply (J7J) by W a (a) on the left. We have 

W Q (a) P Q W(l/a+) W{l/aJ) P a = W a (a) W a (l/a+) W(oZ l ) P Q 
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= [W a {a„)-Q a H{a) H{l/a + ) Q a ) W{l/aJ) P a = I-Q a H{a) H(l/a + ) P a W(l/al) Q. 

= I-Q a H{a) if (l/o+) W(l/al) Q a + Q a H(a) H(l/a + ) (/ - P a ) W^' 1 ) Q a . 
The first two terms combine as 

I-Q a H{a)H{l/{a + aJ))Q a . 

Similarly, 

-W a (a) Q a H(l/al) H(l/a+) Q a = —Q a W(a) P a H(l/al) H(1/<t+) Q a 

= -Qa W(a) H(l/al) H(l/a+) Q a + Q a W(a) (I - P a ) H{l/al) H(l/a + ) Q a . 
The first term here is 

Q a H(a) W(l/*-) H(l/a+) Q a = Q a H(a) H(l/(a- a+)) Q a , 

cancelling the same term above. It follows that W a (a) times equals / plus 

Q a H(a) H(l/a + ) (I - P a ) W(l/al) Q a + Q* W(a) (I - P a ) H(l/al) H(l/<r+) Q a 5 

For convenience we set 

c_ = r _(p), c+ = T+(-p). 

The Hankel operator H(l/a + ) on the left above is the sum of two terms. One is 
c^ 1 e p ® e p while the other is H(u + ). Consider the contribution of the latter. It has 
kernel 0(e~ s(x+y) ), so H(u+) (I - P«) has kernel 0(e~ Sa e~ 5{x+v) ). Since the kernels 
of both H{a) and (J — P a ) W(l/al) Q a are bounded it follows that the contribution 
of the H(u+) summand of H(\/a + ) in the first term is 0(e~ Sa ). Similarly so is the 
contribution of the H(ul) summand of H(l/al) in the second term. All this is 
uniform in p. Therefore with this error the above may be replaced by 

— Q a H(a)e p ® Q a W(l/cr-) (/ - P«) e p + —Q Q W{a){I - P a ) e p ® Q a H(l/a + ) e p . 

C_j_ C— 

Using the general fact Q a W(v) (I — P a ) e p = e tap P a H(v) e p , which comes from 

/■oo roo 

/ V(a-x-y)e ipy dy= / V(-x - y) e ip{a+v) dy, 

Ja JO 

5 In the regular case with the stated conditions the Hankel operators are all compact, so this 
operator is o(l). Similarly for the one obtained by multiplying by W a (cr) on the right. This gives 
J2J almost immediately. With the stronger assumption on a this operator is 0(e~ Sa ). 
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we see that the above equals 



e 



lap piap 



Q a H(a)e p <g> P a H{l/aJ) e p + P a H(a) e p <g> Q a H(l/a+) e p 

giap giap 

= Qa H(a)e p ® P a rj-\ P a H(a) e p (g) Q a C, 

C— 

in our notation. 

This suggests that we add to our first guess a linear combination of tensor products 
with right factors P a i] and Q a (- By symmetry (or taking transposes of the above 
applied to a) the left factors should be P a ( and Q a rj. So it should be of the form 

{aP a ( + bQ a r]) ®r]+ {cP a ( + dQ a r]) <g> Q a (. 

To see if this works, we compute W a (a) P a ( and W a (a) Q a rj. We have 

W a (a) P a ( = P a W(a) P a H(l/a + ) e p = P a W(a) H(l/a+) e p -P a W(a) (I-P a ) H(a^_ 1 ) e p 

= -P a H(a) W{l/<T + ) e p + -^—P a W(a) (/ - P a ) e p (9) 

Zip C_|_ 

1 _ e iap 

P a H(a)e p + — Q a H(a)e p . 



Similarly 



a+(p) ' ' ' 2ipc + 
W a (a) Q a r] = Q a W a (a)r] 



D tap 



Q a H(a) e p + — P a H(a) e p . 



This will cancel the previous extra terms when we add, if our coefficients are such 
that in aP a ( + bQ a r] the coefficient of P a H(a)e p is zero while the coefficient of 
Q a H{a) e p is —e tap /c + , and in c P a ( + dQ a rj the coefficient of Q a H(a) e p is zero 
while that of P a H(a) e p is — e iap /c„. We see that this does occur when 

_ 2ipa_(-p) <T + (p) e 2iap _ (2ipf r_(-p) r+(p) e 2iap 

(X — 



b = 



d = 



A 


A 


(2ip) 2 c_ cr_(-p) e iap 


(2ip) 3 r_(p) r_(-p) e iap 


A 


A 


{2ipfc + a+(p) e iap 


(2ip) 3 r + (p)r + (-p) e iap 


A 


A 


2ipa^(—p) <J + {p) e 2iap 


(2ip) 3 r_(-p) r+(p) e 2iap 
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where 



A = {2ipf c_ c+ - e 2iap cr-{-p) a+(p) 



(2ip) 3 e 



.4 



Thus 



a = r_(-p) r+(p) e iap A, 6 = -r_ (p) r_(-p) A, 
c = _ r+ (p) T+ (-p) A, d— T -(-p) r+(p) e lQp A. 



This tells us that when we define T a by (JHJ) then the first statement of the lemma 
holds. (The reason for the extra factor A in the statement is that the error inherent in 
© gets multiplied by a factor proportional to A.) The second is obtained by applying 
the first to o and taking transposes. This completes the proof of the Lemma. 

Proof of Theorem 1. The theorem follows from the relations of the lemma by using 
the easily checked fact that T a differs from I by an integral operator whose kernel is 



Since the statement of the lemma to Theorem 1 holds uniformly in p, we need only 
compute the p — > limit of T a to obtain its analogue when p = 0. The computations 
simplify if we take r±(0) = 1, which we may do since the theorem for the general case 
follows by applying the special case to the symbol c(£)/t(0). We easily compute that 



For the rest, a tedious but straightforward computation 6 shows that all terms of 
(jHJ) except for the first combine as —B^ 1 times 



[P a H(u + )l-Q a H(^)l-x+ir^(0)}®[P a H(iJl)l-Q a H(u + )l-x-iT'_(0)} + O(p). 



0(a). 



2. Proof of Theorem 2 



.4 



(a + t(r'_(0)-r' + (0))r' + O(p) 



B- l + 0(p). 



(All terms involving p 



cancel.) Now 



W(l/a + )l 



x + W(u+)1 



while H[u + )1 has Fourier transform 




so 



H(u+)1 = -W{u+)l+u + (0) 



W(u+)l-ir' + (0). 



The tedious computation was done by Maple. 
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(Observe that r±(£) = 1 =F i£ w±(£).) Therefore 

W{l/a+)l = x- H(u+)1 - ir' + {0), 

and similarly for the right side of the tensor product. Therefore the tensor product 
equals 

[P a W(l/*+)l + Q a H{ul)l] ® [P a W{l/al)\ + Q a H{u+)1}. 
Thus if we set 

T° = P a W(l/a+) W(l/a-) P a - Q a H(ul) H(u+) Q a 

B- 1 [Q a H(ul)l + P a W(l/a+)l) ® [Q a H(u + )1 + P a ^(1/^1)1], (10) 
then we have 

W a (a)T° = I + 0(Ae~ 5a ), T® W a {a) = I + 0(Ae~ Sa ). (11) 

This gives the statement of Theorem 2 when we check that T° differs from I by an 
integral operator whose kernel is 0(a). 

3. Proof of Theorem 3 

Suppose we want to prove © for a symbol <j\. If we take any other <7o of the 
same form we can write (J\ = <J e p , and we define the family a = a\ = <jq e Xp where 
A G [0, 1]. We assume first that A satisfies the bound in the statement of Theorem 1 
for all A G [0, 1], so that we can use its conclusion and the relation 

-^log det W a (a) = tr Waiay'Waiap). 
dX 

The operator on the right is 

Wa(a)- 1 [W a (a) W a (p) + P a H(a) H(p) P a + Q a H(a) H(p) Q a \ 

= W a (p) + Wa(a)- 1 [P a H{&) H(p) P a + Qa H(d) H(p) Q a \. (12) 
We look first at 

tr Wa(a)- 1 P a H (a) H(p)P a , 

and in particular to the contribution to this of bracketed expression in (JHJ) with its 
factor A. 

First, the factors P a may be removed with error 0(e~ Sa ). Then H(p) may be 
moved around to the left, so the trace of the product = A times that of 

e iap r + (p) r_(-p) H(p) ( ® H(a) V - r + {p) r + (-p) H(p) ( ® H{a) Q a ( 
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-r_(p) r_(-p) if(p) Q a r? ® //(a) ^ + e** r+(p) r_(-p) if (p) QaC ® # (<r) Q^. 

The function if (a) Q a ( contains as a summand H(a) Q a H(u + ) e p (the rest comes 
from the e p <g> e p part oh ff(l/cr + )), and this is 0(e~ Sa ) since H(a) and H{u+) have 
exponentially decaying kernels. Similarly for rj. Therefore with this error we may 
make the replacements 

e iap e iap 

2ipc + p 2ipr + (—p) pi 

e iap e iap 
Q a Tj — > 6_ = 6 

2ipc_ p 2ipr + (p) p 
Next, for any plus function v + the Fourier transform of H{y + )e p is 

I v + \ = ~ v + 

\0 + i(Z-p)) + 0-i(S- p y 

Therefore if we denote the Fourier transform by JF we have in particular 

_ ^(p)- 1 -^- 1 (2ipr_(-p))- 1 + ^- 1 

77 O-i(^-p) O-i(^-p) 

and similarly 

TC = v+iPy 1 -^ 1 = (2tpr + (p))- 1 + a- 1 

0-i^-p) o-KZ-p) 

If we make these substitutions we find that we are left with four tensor prod- 
ucts, only two of which have nonzero coefficients. We show them below with their 
coefficients. 

H{p)e_ p ®H{o)e_ p , f°\ A~\ (13) 

2ipT + (-p)T-{p) 



Next we consider P a H(a) H{p) P a left-multiplied by the Q a if (l/cx_) H(l/a+) Q a 
part of W a (a)~ l . Again, because of the Q a we may make the replacements 

Q a if (1/51) -> el 1 e iap e_ p <g> ep, if (l/o+) Q a -> c; 1 e ia % <g> e_ p , 

the product of these being 



e 2iap 



2ip c_ c + 

11 



6 — p 



and after bringing the H(p) around, the corresponding tensor product is H(p) e_ p 

£ 2iap 



H(a) e_ p with the factor 



2ipT + (-p)r-(p)' 

This cancels the term in (|13|) when we multiply it by A and subtract this one. 

To compute the trace of the term in (fTf|) we note that the factor on the right has 



Fourier transform 



o 1 \ (J i 



since p) = 0. Therefore the inner product of the two factors is 

J_ f V+ 1 ^+ d$ \_ f P_ d , 

27rJ p o-i(e-p) o-z(e+p) * ' 2nJ ( -i(e-p))(o-z(e+p)) 5 

and so the trace with its factor (and the factor A) is 

Ae'°>T +{r) T4-v) P - {P) - p - { - p) . (15) 

hp 

As for the P a W(l/a + ) W(l/a-) P a part of the inverse, this is = the trace of 

W(l/a + ) W(l/a-) H(a) H(p) = W{l/a + ) H(a + ) H{p) 

= -H(l/a+) W(cT + ) H(p) = -H{l/a+) H(pa^). 
The trace (including the minus sign) equals 

±f<p, + Y,?it=±fr-Oo e „ + M, (i6) 

since / p' d£ = 0. 

So (|T3 |) + (|T^j) is the contribution of the P a H(a) H(p) P a term in (JT2j) to the trace. 
For the contribution of the Q a H(a) H(p) Q a term, observe that 

W a (a) 1 Q a H(a) H(p) Q a = Q a W a (a) ~ x H(a) H(p) Q a , 

which has the same trace as 

W a (a)- 1 H(a)H(p)P a . 

Hence we need only replace a by a everywhere in what we derivated above. Then we 
add two and we get 
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2ip 



+ 



2?ri 

And finally we have to add to this 



X -\\p- (log<7 + )'-p + (lo g( 7_)V£- 



tr W a (p) = ^ J p(0 df. 



What is this the derivative of? Since (d/dX) loga = p, for the last term the answer 

is 

Here we use our chosen factorization of a to define the logarithm. For the next to 
last terms, since (d/dX) log<x± = p± we have 

J- J (l0g<7+)' lo g( T_ rfe = y [/>'+ lo g( 7_ + (l g<7+)V] d£, 

so the middle term is d/dX of 

i/ (1 ° g(T+)/ l0ga -^ 
For the first term we compute that since dr±/dX = p±r±, 

d 1 . d . 1 
log A = — log ■ 



dX dX r_(p) r + (-p) e~ ia P - r + (p) r_(-p) e ia P 

= -jp-ip) + P+(-p)) t -(p) t +(-p) e~ iap + (p+ip) + P-(-p)) r +(p) t -(-p) e iap 
T- (p) r+(-p) e~ ia P - t + (p) T-(-p) e ia P 

In the numerator we add (p-(p) + P+(— p)) T +(p) r -(~ P) e * ap to the first summand 
and subtract it from the first. Then the quotient becomes 

-p-(l>) - M-P) -Ae>°°r +{P ) r.(-p) P-W - P-(-P) + P+(-P> - PM , 

lip 

and so the first term is the logarithmic derivative of 

(r + (-p)r.(p)A)-\ 

which is the same as the logarithmic derivative of 

r_(p) r+(— p) e~' tap — T + (p) r_(— p) e iap 
t+(-p)t-(p) 
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Putting this together gives the logarithmic derivative of det W Q (a) with error 
0(e~~ 5a ). Integrating and exponentiating gives 



det W a (a) 
G(<r)° 



C exp j— J (logo+)' loga_ d^J 



r-jp) r+j-p) e~ iQ P - r + (p) r_(-p) 

r + (-p)r_(p) ' (17) 

where C is a constant independent of A. 

This was derived under the assumption that, for all A, A = 0(e Sa ) for some small 
S, which means that 2iap is not exponentially close to any of the values of 

r_(p)r + (-p) 
log^^ : r, 18 

in particular if 2iap is bounded away from these values. Assume that this is so for 
our original a±, and suppose we had taken 

MO 



with gi, g 2 in the right half-plane. The corresponding quotient in (fTSj) is 

gi ~ ?2 - ip 

qi + ip q2 + ip 

Each factor maps the right half-plane to the lower half-plane, so with proper choice 
of qi the ratio can be made arbitrarily close to the ratio in (fT5j) coming from <j\. The 
values of (|18|) for o\ depend linearly on A, so for this choice 2iap will be bounded 
away from the values of (J18)) for all A G [0, 1]. Therefore (|17j) holds for some C. We 
have p, Prop. 5.2] 



det W a (a ) = e 



-a(g 1 +q 2 ) 



(qi + ip) (q 2 + ip) e tap - {q x - ip) (q 2 - ip) e 
2ip (qx + q 2 ) 



-a(pi+p2)/2 



and a computation shows that for this o"o (fTTj) holds when C = e iap . 

This proves the theorem for our symbol <j\ under the condition that 2iap is 
bounded away from the values of (|18|). To remove the condition, suppose that it 
is very close to one of these values. Replace 0\ by o\ e Xp , with p now any symbol such 
that 

P-(p) P+(~P) 



P+(P) P-(-p) 



7^1, 
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and let A run over a little (but not too little) circle around zero. Then the corrre- 
spoinding values of (fTSj) will run over circles about the original ones, and so 2iap will 
be bounded away from them, uniformly in A. By analyticity of the two sides in our 
formula, their values at A = are obtained by integrating with respect to A over the 
circle. We deduce that the difference of the two sides, which is exponentially small 
uniformly for A on the circle, is also exponentially small for A = 0. 



4. Proof of Theorem 4 



One can check that the right side of © is the limit of the right side of (0) as 
p — ► 0. Unfortunately it is not clear from our derivation that (JHJ holds uniformly in 
p, although it surely does. Fortunately the derivation of Theorem 4 from Theorem 2 
is a lot easier than the derivation of Theorem 3 from Theorem 1, so we proceed with 
it. 

We start in the same way, with (|12|). As before, the first term has trace 

and we proceed to the contribution of 

trW a (a)- 1 P a H(a)H(p)P a . 
The contribution of the P a Wi^aZ 1 ) W(aZ 1 ) P a part of the inverse is 

JL/^(W««. us) 

as before. The contribution of the rest of W^a)' 1 is the trace of 

-Q a H(ul) H{u+) Qa-B- 1 [Q a H(€l)l+P a W(al l )l]®\Q a H(u + )1+P a W^ 1 )!} 

right-multiplied by P a H(a) H(p) P a . In computing the trace we may, as before, 
move the H(p) P a around to the left. As we saw earlier, any trace involving two 
Hankel operator with nonsingular symbol with the operator Q a between them will 
be 0(e~ Sa ). With similar error we may remove the factors P a . So with this error the 
trace reduces to that of 

-B- 1 H(p) W^al 1 )! ® H{a) W^- 1 )!. 

The Fourier transform of the second factor is 



a 



°- 1 \ a + 
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while the Fourier transform of the first factor is 

The p may be replaced by p+, and the result after replacing £ by — £ is 



P- 



al 1 



0-if 



Notice that the expression in parenthesis has the singularity (0 — i£)~ 2 - When we 
multiply by the preceding and integrate we may remove the minus subscript, and we 
see that the trace in question equals 



So the contribution of the P a H(a)H(p)P a factor is — iB~ x p'_(0), which is to be 
added to ([19]). Similarly (or by applying the preceding to a) we compute the contri- 
bution of the Q a H{&) H(p) Q a term. Adding the two gives 

(log<7 + )'d£ - ^ fp+ Qog<7_)'d£ + iB- 1 (p'_(0) - p' + (0)). 

Thus the logarithmic derivative of the determinant equals 

7^ / P(0 rf£ + — / p- (log cr + )'di+^-. I p' + log a_ d£ 

(p'_(0) - p' + (0)) + 0(e- Sa ). 

The first term is the logarithmic derivative of G(a) a , and the next two terms combine 
as the logarithmic derivative of 

6XP {dri / ( loga+ y loga - ' 
as we saw in the last section. Since (see footnote H]) 

B = a-i (\og^ + \{p + -p^ (0), 

the last term above is the logarithmic derivative of B. This proves the theorem with 
a constant factor C on the right side. Taking the special case er (£) = £ 2 /(l + £ 2 ) 
shows that C = 1. 
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5. Comparison with the formula of [lj 

We show here that (jHJ) agrees with the result of jl.. Take p > 0. We define 
two contours, both going from — oo to oo. Contour C\ has an indentation above — p 
and below p, while C 2 has an indentation below — p and above p. Each has its own 
geometric mean Gj(er) and constant Ei(a). These are defined by integrals over 
with integrands involving the associated factorizations. Our factorization a = o_ a + 
is associated with C\. Thus G{p) = G±(a), E(a) = E\(cr), and if we think of the right 
side of (jSJ as a sum then G(a) a times the first summand is exactly Gi(a) a E\(a). 

For the other terms we first express our exponential factor Ei(a) in terms of E 2 (a). 
We have 

log£i(a) = ^-J Ci (log((e + p)r + ))'(log((e-p)r_))de 

1 roo-iS . s 1 . s 

= 2" /_ oo _. (5 ( M(f + P) ( M(f - p)t-)) d£ - log(-2p r_(-p)), 
log£ 2 (a) = ^/ ca (log((e-p)r + ))'(log((e + p)r_))de 
= — / (log £-p r+ ) (log £ + p r_ W-k)g 2pr_ (p , 

Z7TZ J-oo-iS v 7 v ' 

since for the first we pass the pole at £ = —p and for the second the pole at £ = —p. 
For the two resulting integrals, both over the contour in the lower half-plane, the first 
one minus the second equals 

-L /(logr+yiogi^de+r^ / (log log((£ -p)t-K 

Z7TZ J £+P 27TZ j \ c, — P J 

\ogr + (log d^ + ^- ( (log log((£ -p)r_)de 



27TZ J V ^ + P 

= - log ^+ (p) + log r+ ( -p) . 
Putting these together shows that 

Ei{a) ■; -—h 2 (a). 

r-{-p) r + {p) 

Finally, we compute that G 1 {a)/G 2 {a) = e~ 2ip . Thus e 2iap G(a) a = G 2 (a) a , so 
G(a) a times the second summand in our formula is exactly G 2 (a) a E 2 (a). Hence our 
result can be stated 

det W a {&) = G x {a) a (Ex(a) + 0(e' Sa )) + G 2 (a) a (E 2 (a) + 0(e~ Sa )), 

which, except for the size of the error terms, agrees with pQ. 
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6. Another proof of Theorem 2 

As before, we may assume r_(0) = t + (0). We begin by multiplying W a (a) on the 
right by T°, given by |T0|) . For the product with P a W(l/a + ) W(l/a-) P a we obtain 
as before 

W a (a) P a W(l/v+) W(1/*J) P a = I-Q a H{a) H(l/a+) P a W(l/vl) Q a , 
but now we write 

H(a) H(l/a+) P a W{l/al) = H(a) (1 ® 1) P a W(l/al) + H(a) H(u+) P a W^ 1 ) 

= H{&) (1 ® 1) P a W(l/al) + H(a) H(u+) W(l/al) 

= H(d) (1 <g> 1) P a W(l/al) + H(a) H(u + /aJ). (20) 
This is to be left- and right- multiplied by Q a . 

Next, for the product with Q a H(ul) H(u+) Q a we heve 

W a (a) Q a H{ul) H{u + ) Q a = Q a W{a) P a H{u_) H{u + ) Q a = Q a W{a) H{u_) H{u + ) Q c 

and we compute 

W(a) H{ul) H(u+) = [H{aul)-H{a) W{u-)} H(u+) = H(aul) H(u+)-H(a) H(u- u+ 
= h(<t(±- ) H(u + ) - H{d) H (u + (— 1 ^ 



= H (oil?) " + H{d) " (orie) - H{d) HMa - ] - 

Here we used the fact that cr(0) = 0. The last term here cancels the same term 
appearing in (|20p . The first two terms we may write as 

Since 

— — + — — = 2tt5 = 1, 
- if - if 

these terms combine as 

H(a)(l®l)H(u+). 
Putting together what we have done gives 

W a (a) [P a W(l/a+) W{l/aJ) P a - Q a H(l/al) H(l/a+) Q a \ 
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= I-Q a [H(a) (1 ® 1) P a W(l/al) + H(a) (1 ® 1) if (u+)] Q a . 
For the right side of the first tensor product we use Q a l = P a l to see that 

(P a W{l/al) Q a )'l = Q a W(l/a_) P a l = P a W(l/al) Q a 1 

= P a W{l/al) P a l = P a W{l/al) 1. 
Thus the right side above equals 

I - Q a H(a)l ® [P a W(1/5I)1 + Q a H(u + )1}. 

Next, we will not make a guess as we did in the proof of the lemma since we really 
know the answer. So let us just continue by computing 

W a (a) [Q a H{ul)\ + P a W{a- l )l] = Q a W{a) H{u_)l + P a W{a) P a W{cj- 1 ) P a 1. 

The first summand equals 

Qa H {^{J~_- qT^) ) # (*) = Q Q if l-U-(O) Q a if (cr) 1 

The second summand equals 

P a W(aJ) P a l -aQ a if (or) ff (1/(7+) Q a . 

Using if (l/cr + ) = 1 <g) 1 + H(u+) and 1 ® P Q 1 = a 1 ® 1 we see that this is = to 

P a W(a.) P a l -aQ a H(a)l - Q a H(a) H(u+)1. 

Since if (u + )l = — W(u + )1 + u+(0) when we combine the terms we get 

P a W(a_) P a l -aQ a H(a)l + Q a H(a) W^ 1 ) 1 - ( u+ (0) + «_(0)) Q a H(a)l. 

Since, as already observed, u±(0) = =Rt±(0), the coefficient in the last term equals 
—% (r+/r_)'(0) = i (o+/er_)'(0). Therefore the above equals 

P a W>_) P a l + Q Q ff (a) ^(a; 1 ) 1 - ff Q Q ff ((x)l. 

Let us show that the first two terms cancel. We have if (a) W{a^_ r ) 1 = H(al) 1. Since 
er_(0) = this has Fourier transform <rl/i£, which is minus the Fourier transform of 
W(al) 1. Therefore the second term above equals 

-Q a W(al) 1 = -Q a W(al) P a 1 = -P a W{aJ) Q a 1 = -P a W>_) P a 1. 
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We computed above what we obtain when we left-multiply the first two terms in 
(JTUJ) by W a (a). What we just computed shows that if we do the same with the last 
term and add, what results is / + 0(e~ Sa ). Similarly for multiplication on the right. 
(Or apply the preceding to a and take the transpose.) This gives (fTTj) and hence the 
statement of the theorem.. 
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